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Abstract 

Using a mean-field BCS-like approach on the bidimensional extended Hubbard 
Hamiltonian we calculate the superconducting transition temperature T c as a func- 
tion of the hole content n/j. This method can be used to determine the critical 
temperatures T c either for the extended — s wave or d wave symmetry. We can also 
describe the pressure effects on T c by introducing the assumption that it induces a 
change in the magnitude V of the attractive potential. This assumption yields an 
explanation for the intrinsic term, and together with the well known change in rih, 
we set the critical temperature as T c = T c (nh(P), V(P)). Then, we obtain a general 
expansion of T c in terms of the pressure P and the hole content n^. We apply this 
expansion to the Tlo^Pbo.5Sr2Cai- x Y x Cu20'r system (0.0 < x < 0.35) and our 
results provide a good fitting for the experimental data for the variation of T c with 
the pressure P and n^. 
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1 Introduction 



In spite of the advance attained in the research of the high temperature su- 
perconductors (HTSC), a great number of questions related to them remain 
to be answered. In order to obtain a better understanding of these new ma- 
terials many experiments were performed, but the different experiments were 
interpreted by different theories. One of the most fundamental evidence that 
these experiments showed was the interplay among critical temperature T c , 
hole content rih, and pressure P, since one of the effects of the application of 
external pressure is a change in the hole content in the Cu02 planes [1], 
which results in a large change in T c . The large increase of T c by the appli- 
cation of pressure [2] has been the motivation for a vast number of works on 
the study of the effects of pressure on the critical temperature T c of the high 
temperature superconductors [3-6], not only to achieve higher T c but also to 
obtain a hint for chemical substitution which originates new materials and 
yields a better understanding of their physical properties. 

Thus, in the detailed study of the dependence of the critical temperature 
T c of the HTSC on the pressure P, it has been suggested that the derivate 
dT c /dP may be strongly influenced by a pressure-induced alteration in the 
density of hole-like charge carriers rih on the CuC>2 planes [7]. This proposal 
was also supported by experimental observations of Hall constant's pressure 
dependence [1] that confirmed the existence of a growth in rih in the Cu0 2 
planes when the pressure is applied [8]. Therefore, there is a pressure induced 
charge transfer (PICT) from the block layers into the CuO<i planes, that is, 
n-h = rih(P)- Since one expects dT c /dP = at the optimum doping, and 
T c (n op ) generally increases with P, it is well accepted that the pressure also 
increases T c by a different mechanism of "intrinsic" origin [3,5,9], and such 
intrinsic term is a characteristic of each family of compounds. Therefore, we 
can separate these two effects by considering [3,5,9,10]: 

dTc = dri m dun m 

dP dP dn h dP' U 

Here the first term, dT^/dP, represents the intrinsic contribution, not related 
to the charge transfer induced by pressure drih/dP. Thus, the change of T c 
with the hole content is given by dT c /drih- 

In this paper we consider the intrinsic contribution as originated from a vari- 
ation of the magnitude V, of the attractive potential that forms the Cooper 
pairs in the superconducting phase with the applied pressure P, in such a 
way that V = V{P). Similar proposition was previously made by Angilella 
et al. [9] and de Mello et al [3,5]. As discussed in Ref. [3,4], as V increases 
with the applied pressure, the zero temperature superconductor gap also in- 
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creases. This hypothesis has been indirectly verified in the experimental study 
of polycristallines samples of HgQ S 2Reo,i 8 Ba 2 Ca 2 Cu 3 Os+8 [11]. 



2 The Method 

In order to study the dynamics of the hole- type carriers in the superconducting 
phase as well as the normal phase with correlations and the basic attractive in- 
teraction, we consider a two dimension extended Hubbard Hamiltonian [3,5,9] 

H=- tijcl a C ja + Uj2 n H n H+ J2 V ij C ia c ]a' C jT' c ia, (2) 

<C«i>o" i <ij>aa' 

where tij is the hopping integral between the nearest-neighbour and next- 
nearest-neighbour sites i and j; U is the on-site correlated repulsion and Vy 
describes the attractive interaction between nearest-neighbour sites % and j. 

In the superconducting phase the paired states form a condensate separated 
from the single-particle states by a "gap" (order parameter). Using a BCS- 
type mean-field approximation [12] the self-consistent gap equation, at finite 
temperatures, is given by 

^"X^tanh^, (3) 

with 

£k = y/el + Al (4) 

Here, Vkk 7 is the interaction potential which contains the repulsive and attrac- 
tive potential, and may be written in a separable form [9,13] 

Vkk' = U + 2V cos(k x a) cos(k' x a) + 2V cos(k y a) cos(A;^a), (5) 

where it was considered a square lattice of lattice parameter a and equal 
coupling constants V along both directions in the Cu-0 planes, in the context 
of mean field approximation. 

A tight-binding approximation may be employed to dispersion relation in Eq.4, 
which yields 

£k = — 2t(cos(k x a) + cos(k y a)) + 4t xy cos(k x a) cos(k y a) — /i. (6) 
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This equation has been previously considered [3,5,9,13], but here it was consid- 
ered identical hopping integrals along both directions in the Cu-0 planes for 
the nearest-neighbour (t x — t y — t) and a different one for the next-nearest- 
neighbour (txy); n is the chemical potential. 

Also, using the same BCS-type mean-field approximation applied to the gap 
(Eq.3), we obtain the hole-content equation [14] 



where < rih < 1. 

Following the steps of Ref [9], one observes that the substitution of the poten- 
tial (Eq.5) in the gap equation (Eq.3) leads to the appearance of a gap with 
two different symmetries 



where the plus sign is for extended — s wave and the minus sign, for d wave 
symmetry. The anisotropy of the gap is one of the characteristic distinction 
between the usual low temperatures superconductors and the HTSC. Many 
experiments and theoretical calculations suggest that the HTSC exhibit a pair- 
ing symmetry different than the usual superconductors [15], and in some cases 
even a mixture of different symmetries must be considered [16]. Thus, the de- 
termination of the symmetry of the gap is the first step in the identification of 
the pairing mechanism of the charge carriers and the subsequent development 
of a microscopic theory for the HTSC. 

To obtain the phase diagram, the gap equations are solved numerically, in the 
limit T — > T c , together and self-consistently with the density of hole carriers 
n h (Eq.7). In these calculations it was considered that there were no symmetry 
mixture near the critical point [9,16]. 

The values of the hopping integrals t and t xy used in the calculation of the 
phase diagram of the Tl 5 Pbo^Sr 2 Cai^ x Y x Cu20i series are taken from A.R.P. 
E.S. measurements for the Bi2212 [9], since it is a HTSC with similar structure, 
having two Cu02 identical layers per unit cell [17], although this does not 
guarantee that the real t and t xy have this values and that is why we have 
tried a small variation around the A.R.P.E.S results. The chemical potential fx 
is calculated self-consistently. The values of U and V are considered adjustable 
parameters, in order to reproduce the T c x rih curve close to the experimental 
data. Our method to derive the phase diagrams is similar to the approach 
based on a change in the zero temperature gap [18]. On Fig. 1 we show the 
phase diagram with the theoretical curves from the numerical calculation for 





A k (/i,T) 



A max (n,T)[cos(k x a) ± cos(fc„a)], 



(8) 



4 



the extended — s and d wave symmetries, together with the experimental data 
for the Tl 5 Pb , 5 Sr 2 Cai_ x Y x Cu2O- 7 [19]. 



90.0 



T (K) 



40.0 





Phase Diagram 




U = 0.052eV t = 0.05eV 


1 


1 i 1 i 

▲ ref 

V = -0.0875eV t xy =0.023eV 

V = -0.057eV t xy =0.0275eV 


A" 




r 


d-wave 


// A 




|i extended-s 


\ '; i , i 



0.00 



0.20 



0.40 



0.60 



0.80 



1.00 



Fig. 1. Phase diagram showing the numerical calculation 
for the extended — s and d wave symmetry together with the experimental data 
of the Tlo^Pbo.sS^Cai-xYxC^Of series taken from Ref. [19]. 

Fig. 2 illustrates the behavior of the phase diagram, for the extended — s 
symmetry, when V is varied. Both symmetries exhibit a similar dependence 
on the potential V. 



Notice that although there is a strong tendency toward a rf-wave symmetry 
on many experiments the HTSC, our results, with the parameters t and t xy 
fixed, show that the extended — s symmetry seems to reproduce well the phase 
diagram of the Tl 05 Pb j J Sr 2 Ca 1 _ x Y x Cu2Oi series in a better agreement with 
the experimental points than the d wave one. This is in accordance with the 
results of Schneider et al. [13]. 
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Fig. 2. Phase diagram of the TlQ^PbQ^Sr2Ca\- x Y x Cu20y series showing the effect 
of a change in the magnitude V of the attractive potential. 

Thus, the study of the phase diagram yields the best values £=0.05eV, t xy —0.023eV, 
[/=0.052eV and V=-0.0875eV for the extended - s symmetry and t=0.05eV, 
t xy =0.0275eV, £/=0.052eV and \/=-0.057eV for the d wave symmetry which 
will be used in the next section. 



3 The Pressure Effects 



On Fig. 2 we show that T c is quite sensible with respect to a change in the 
magnitude V of the attractive potential. Therefore, as it was previous men- 
tioned, we suppose that the pressure induces a change in the attractive po- 
tential, in such a way that V = V(P). In the absence of informations of how 
V effective behaves with respect to the pressure P, we assume a linear depen- 



6 



dence [3,5,20], which should be the case since the pressure causes very small 
structural changes in the lattice constants [9]. Therefore 



V(P) =V + AV(P). (9) 

In the same way, for the density of carriers, we have 

n h {P) =n h + An h (P). (10) 



Here, V and n h are the P = values of V(P) and n h (P). The pressure 
dependence is defined as 

AV(P)= Cl P, An h (P)=c 2 P (11) 



where C\ = (dV/dP) and c 2 = (dn h /dP) are parameters usually determined 
from the experimental data. 

Our previous discussion leads to T c (rih,P) = T c (nt l (P),V(P)). Thus, we can 
estimate T c for a compound with a certain value n^, and under pressure P, 
using an expansion of T c {n h , P) in terms of P. Consequently, we have 



where 



T c (n h , 0) is the critical temperature for P—0. The Eq. (12) can be written as 

T C K,P) = J> 2 ^, (14) 



with 



d 



d 



az = { Cl dv + C2 d^ h 



T c (n h (P),V(P)). 



(15) 
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One can see by considering the case z — 1 in the above equation, and com- 
paring with Eq.(l), that 

dJi = dvdn ( . 

dP dPOV 1 ' 



which clearly shows that, in our model, the intrinsic contribution is related to 
the variation of the magnitude V of the attractive potential with the pressure 
P. Now, the analytical expressions for the coefficients a z will be derived below 
as a function of the parameters c\ and c 2 , for the values of rih and P. 

The first coefficient (z — 1) for the above expression is given by 

/ dT c dT c \ 

ai = [ Cl dv +C2 d^ h l (17) 



Restricting ourselves to small changes in V, we may approximate 



dT c AT C 

W^AV' (18) 



where the horizontal bar denotes a "mean" over the values of T c obtained 
direct from the phase diagram as a function of V, as it is shown in Fig. 2. 
This "mean" is realized for each value of n. To estimate dT c /dnh we may 
either use an experimental value, or use directly the curves of Fig. 2 or a 
phenomenological universal parabolic curve [19]: 

Tc = T max [i-A(3{ nop -n h fl (19) 



and, therefore, an analytical expression for dT c /drih is obtained. Thus, we 
derive the following expression for the coefficient a.\ 



= c iTT7 + c 2 2(4(3)Tr x (n op - n h ), (20) 



and CK2, which is given by 

( d 8 \( 8T dT c \ 

a2= { Cl W + C2 ^j{ Cl dv +C2 d^ h )' (21) 



and results in 



ATmai 

a 2 = 2(4P)^ r (n op - n h ) Cl c 2 - 2(4/3)T c m ^, (22) 
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where ( AT ^J X ) is an empirical mean derived numerically by inspection on Fig. 
2 at the optimum doping compound n op . Following along the same procedure, 
the third coefficient becomes 



■max 



ATra 

«3 = -4(4/3)^^4 (23) 



As one can notice the coefficients a z , for z > 4, are all nulls. 



4 Comparison with experimental data 



In order to study the effects of the pressure on the critical temperature of 
any given family of compounds, we need first to calculate the constants C\ 
and c 2 , which yield the coefficients a z . Thus, we can see that the constants 
Ci and c 2 are determined by how V and n^ changes with the pressure (Eq 
(11)). Hall effects can give directly drih/dP as in the case of Bi2212 [8] 
and therefore c 2 [4-6]. Since these measurements are not available for the 
Tl 5 Pb 5 Sr 2 Cai_ x Y x Cu 2 7 , we can obtain C\ and c 2 by comparing with the 
curves of T c x P for two values of n^. In fact Wijngaarden et al. [17] also esti- 
mate the value of c 2 from their experimental data for the Tl _ 5 Pb _ 5 Sr 2 Cai- x Y x 
Cu 2 7 . 

Since at low pressures only the linear term comes into play, the higher order 
coefficients can be ignored. Thus, ot\ becomes the slope of a linear approxi- 
mation from the initial points of the T c x P curve, for a given rih compound. 
Starting with nh=n op , in spite of not being strictly necessary, we determine c\ 
using the estimated ol\ in Eq.(20), as long as at n op the charge transfer term 
vanishes. To determine c 2 we estimate «i from the T c x P curve of an n^ ^ n op 
compound, and use again Eq.(20). Once these two constants are determined, 
the a z coefficients for any n h value can be calculated (Eq.(15)). 

For the Tl , 5 Pb . 5 Sr 2 Ca 1 ^ x Y x Cu2O 7 [17,19] series we used the experimental 
values for the rih = n op = 0.15 and ^=0.20 compounds to calculate c\ and 
c 2 . The estimated ct\ were 1.5K/GPa and OAlK/GPa for n op and ^=0.20, 
respectively. The phase diagram parameters (^y ) taken from Fig. 1 and 2 
and used in the calculations were 2900fT/ eV for the extended — s wave and 
2300K/eV for the d wave symmetry for the n op compound and 2700-ft'/el / 
for the extended — s wave and 2600iT/el / for the d wave symmetry for the 
7^=0.20 compound. The parameters of the parabolic curve (Eq.19) used were 
(3 = 3.5 and T c ma:E =105.1K. Therefore, the resulting constants were c\ = 5.17 x 
lCT 4 eVy GPa and c 2 = 6.7 x lO^GPa^ 1 for the extended — s wave and cf = 
6.52 x lQ~ A eV/GPa and c d 2 = 8.7 x lO^GPa- 1 for the d wave symmetry. 
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On Fig. 3 and 4 we present our results in comparison with the experimental 
data. One observes that, for the optimum doping, and the compounds near 
the optimum doping, we have a rather good result. For the optimum com- 
pound (nh— 0.15) in particular, we had the better result, for both symmetries, 
with very good agreement with the experimental data [17,19]. For ^=0.075, 
0.20 and 0.25 we can reproduce well the low temperature data and obtain a 
qualitative agreement on the high pressure region. 




40 1 1 1 1 1 1 — 1 

5 10 

P(GPa) 

Fig. 3. Numerical calculation for the extended — s (filled lines) and d wave (dashed 
lines) symmetries, in the overdoped region, together with the pressure experimental 
data of the TlQ^,PbQ^Sr2Ca\- x Y x Cu20-j series taken from Ref. [19]. 

The value of the constant c 2 (= 6.7 x lO^GPa' 1 ), which means the rate of 
charge transfer by the pressure (^p), for the extended — s symmetry is close to 
the experimental results for the YBCO compounds [10], and is smaller than 
the experimental Hall coefficient value for the Bi2212 [8]. The value estimated 
experimentally by Wijngaarden et al. [17,19] for the Tl . 5 Pbo.5Sr2Cai_ x Y x Cu2 
7 series was dn h /dP = 3.9 x 10 _4 GPa _1 . This discrepancy can be at- 
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Fig. 4. Numerical calculation for the extended — s (filled lines) and d wave (dashed 
lines) symmetries, in the underdoped region, together with the pressure experimen- 
tal data of the TlQ^Pbo^Sr2Ca\- x Y x Cii20j series taken from Ref. [19] 



tributed to the parameters which are difficult to be estimated with the 
very few experimental points, what introduces some uncertainties in the cal- 
culated phase diagrams. This discrepancy is also the reason why our curves is 
in good agreement with the experimental data up to lOGPa, but starts to de- 
viate above this value. Another evidence that our main source of error comes 
from dnh/dP is the fact that the optimum compound, with ^=0.15, which is 
the one with drih/dP=0, is in much better agreement with the experimental 
results than the others, with different n^. 
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5 Conclusions 



In conclusion, to relate the pressure P, the density of holes n h and the attrac- 
tive potential V we used in the calculation a BCS-type mean-field approach on 
the bidimensional extended Hubbard Hamiltonian with the potential V as an 
adjustable parameter. Our method has general application, and can therefore 
be used in any compound under pressure. We showed that these calculations 
can describe reasonably well the pressure effects on the critical temperature 
T c for the series Tl 05 Pb 05 Sr 2 Ca 1 ^ x Y x Cu 2 O 7 , for the four measured doping 
values. 

Our results for the extended — s wave symmetry give a better quantitative 
agreement with the experimental data. But the ci-wave calculations are also 
in qualitative agreement with the experimental data, showing that the order 
parameter symmetry is not decisive to study the pressure effects. 

Thus, our method provides a very interesting novel explanation to the intrin- 
sic term, which is related with the change in the potential V, in agreement 
with a previous hypothesis that the pressure increases the zero temperature 
superconducting gap [3,4] and with some experimental results [11]. 
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